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Abstract 



The radiative effects by fiigh-redshift ultraviolet radiation background (UVB) 
fieats tlie gas in the IGM and ehminate the neutral hydrogen and helium that 
dominate the cooling of the primordial gas. We investigate the role of the radiative 
effects for the temporal evolution of the reionization fraction by using cosmological 
Smooth Particle Hydrodynamics (SPH) simulations. 

We find that the increase of photo-ionization and photo-heating rates due to 
optical depth effects results in a significantly contribute to the heating of the IGM 
before and during the reionization. The main effect of the UV radiation spectrum 
on the temporal evolution of the ionization fraction is given by the value of the 
reionization redshift, Zre, and the redshift interval, Az, in which the reionization is 
completed. 

We evaluate the effects of the UV radiation background on the CMB angular 
power spectrum taking into account different temporal evolutions of the ionization 
fraction. We find that for reionization models with degenerated CMB temperature 
anisotropy power spectra, the radiative mechanisms leave distinct signatures on the 
E-mode polarization power spectrum, at large scales (I < 50). We show that through 
E-mode CMB polarization power spectrum measurements, the Planck experiment 
will have the sensitivity to distinguish between different reionization histories even 
when they imply the same optical depth to electron scattering and degenerated Ct 
power spectra. 

This work has been done in the framework of the Planck LFI activities. 
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1 Introduction 



The detailed study of the intergalactic medium (IGM) is fundamentally im- 
portant for understanding the large-scale structure properties and the galaxy 
formation process. 

At the epoch of reionization the collapsed objects began to influence the dif- 
fuse gas in the IGM and rendered it transparent to the ultraviolet photons. In 
the hierarchical models of structure formation, an early UV radiation back- 
ground (UVB) is produced by the first collapsed halos with masses near the 
cosmological Jeans mass. In order to virialize in the potential wells of the 
dark matter halos, the gas in the IGM must have a mass, Mb, greater than 
the Jeans mass, Mj, which is ~ lO^M© at a redshift z ~ 30, corresponding to 
a virial temperature T^j^ ~ lO^K. Photoionization by the high-redshift UVB 
heats the low density gas in the IGM before it falls into the dark matter wells, 
strongly reducing the fraction of neutral hydrogen and helium that dominate 
the cooling of the primordial gas at temperatures of T^ir. 
The primordial gas dynamics in the photoionized IGM gives important in- 
sight into the epoch of the reionization and the end of the "dark ages" of the 
universe, believed to have been caused by the light from the first generation 
of massive stars (see, e.g., the review by Barkana and Loeb (2001) and the 
references therein). 

From the observational point of view, the study of the reionization process it- 
self as well as the properties of the sources driving it is challenged by a variety 
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of observational probes. A powerful observational probe comes from the Lyo; 
absorption spectra of the high redshift quasars (Becker et al., 2001; Fan et 
al., 2002; White et al., 2003; Fan et al., 2004) showing that all known quasars 
with z > 6 have a complete Gunn-Peterson (GP) trough and a rapid evolving 
hydrogen neutral fraction compatible with the final stage of reionization. This 
conclusion is strengthen by the study of the proximity effects around high 
redshift quasars (Wyithe and Loeb, 2004; Mesinger and Haiman, 2004). Also, 
the high value of the IGM temperature inferred from the study of Lya forest 
at 2; ~ 4 indicates that the hydrogen neutral fraction changed substantially at 
z<9 (Theuns et al., 2002; Hui and Haiman, 2003). 

An other powerful observational probe is represented by the high value of 
the electron scattering optical depth, Tes, inferred from the Cosmic Microwave 
Background (CMB) anisotropy measured by the WMAP experiment (Kogut 
et al., 2003; Spergel et al., 2003; Verde et al., 2003) which requires reionization 
to begin at z > 14. 

In order to reconcile these observations a number of semi-analytical models 
have been developed, the most intriguing possibility being the so called "dou- 
ble reionization" for which the globally-average ionization fraction decreases 
with cosmic time over a limited period. The resulting reionization histories 
display a wide range of features extended over a long redshift interval. Some 
of these works (Cen, 2002; Wyithe and Loeb, 2004) agree on the necessity of 
a first generation of metal-free stars (Population III) with heavy mass func- 
tion to produce the large observed value of Tes- Other studies (Ciardi, Ferrara 
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and White, 2003; Somerville and Livio, 2000) found that the metal enriched 
stars (Population II) are able to reionize the universe enough early to produce 
such high values of Tes- These works assume that the total volume fraction 
of the ionized regions were driven by the ionizing sources located in the dark 
matter halos. Typically, these scenarios require a number of assumptions on 
the dynamics of the stellar evolution. Recent works (Scannapieco et al., 2003; 
Furlanetto and Loeb, 2004) show that the transition redshift at which the star 
formation switches between the two modes must occur over a large redshift 
interval and only after the reionization is complete. In this scenarios, the ob- 
servables characterizing the global reionization history of the universe can be 
reconstructed with a large range of input parameters. 

Other physical mechanisms that can explain the extension of reionization over 
a long redshift interval are radiative, the most important being the photo- 
ionization heating. A recent work (Furlanetto and Loeb, 2004) critically exam- 
ined the plausibility of different mechanisms, showing that double reionization 
requires a rapid drop in ionizing emissivity over a single recombination time 
that can be obtained with unusual choices of the physical parameters. 
Cosmological N-body simulations provide a completely different approach to 
model the reionization process (Ostriker and Gnedin, 1996; Gnedin and Os- 
triker, 1997; Gnedin, 2000; Ciardi, Stoehr and White, 2003), offering the ad- 
vantage over the semi- analytical models to be able to fully account for the 
dynamical evolution of the matter contents in the universe. The main limi- 
tation of N-body simulations is related to their resolution that turns to be 
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critical for the reionization studies. At redshifts z > 3 the quasar population 
declines (Madau, Meiksin and Rees, 1999) and it could not be able to reionize 
the gas in the IGM. The bulk of the ionizing photons at high redshifts should 
be then produced by an early population of stars (see, e.g., Haiman, Abel and 
Rees, 2000, Barkana and Loeb, 2000, Cen, 2002). 

In the hierarchical theories of structure formation the low mass objects col- 
lapse first merging then into progressively larger systems. In this scenario, 
small mass objects (M ~ 10^ — IO^Mq) are the main contributors at high 
redshifts while their role becomes much less important at low redshifts in 
comparison with that coming from objects with higher masses (M > lO^M©). 
Thus, numerical simulations are subjected to the difficulty to be able to resolve 
objects in a large dynamic range with enough resolution. 

A number of works have also considered the implications of different reion- 
ization mechanisms for the CMB angular power spectra (see e.g., Bruscoli, 
Ferrara and Scannapieco 2002; Holder et al., 2003; Haiman and Holder, 2003; 
Naselsky and Chiang, 2004). 

Although the electron optical depth to the Thompson scattering of the CMB 
photons, is an important cosmological parameter (see, e.g., Kaplinghat et al., 
2003; Hu and Holder, 2003), the CMB angular power spectrum contains more 
information about the reionization process than the optical depth integrated 
over the whole ionization history. 

In this paper we study the role of the radiative effects for the time evolution 
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of the global ionization fraction at sub-galactic scales (M ~ 5 x lO^M©) by 
using N-body cosmological hydrodynamical simulations. 
We evaluate the effects of the UV radiation background amplitude and shape 
(in both frequency and redshift) on the features of the CMB anisotropy and 
polarization angular power spectra and address their detectability by ongoing 
and future CMB anisotropy experiments. 

The paper is organized as follows: in § 2, we compute the UV radiation back- 
ground flux as solution to the cosmological radiative transfer equation and 
discuss the radiative mechanisms that can cause the global ionization fraction 
to evolve non-monotonically with the cosmic time. We present the various 
reionization histories obtained for the considered UVB flux models in § 3. 
In § 4 we compute the CMB angular power spectra, in both temperature and 
polarization, for the various reionization histories by taking into account the 
radiative mechanisms. The possibility to distinguish among different reioniza- 
tion scenarios with the forthcoming Planck experiment is also discussed in 
this section. We draw our main conclusions in § 4. 

Throughout we assume a background cosmology consistent with the most 
recent cosmological measurements (Spergel et al., 2003) with energy density 
of fljn — 0.27 in matter, flf, = 0.044 in baryons, fl\ = 0.73 in cosmological 
constant, a Hubble constant of Hq—12 km s~^Mpc~^, an rms amplitude of (jg = 
0.84 for mass density fluctuations in a sphere of radius 8h~^Mpc, adiabatic 
initial conditions and a primordial power spectrum with a power-law scalar 
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spectral index Ug — 1- 

2 Building up the UV radiation background 

Reionization is an inhomogeneous process that proceeds in a patchy way. The 
radiation output associated with the collapsed halos gradually builds up a 
cosmic UVB. 

At early times, most of the gas in the IGM is still neutral and the cosmological 
Hn regions around the individual sources do not overlap. At this early stage 
the gas in the IGM is opaque to the ionization photons, causing fluctuations 
of both ionization fraction and UVB intensity from region to region. At the 
reionization redshift, Zre, the Hn regions surrounding the individual sources in 
the IGM overlap. As shown by numerical simulations (see, e.g., Gnedin, 2000, 
Barkana and Loeb, 2000), at Zre the UVB intensity increases sharply and the 
baryonic gas is heated to a temperature Tjgm — lO^K. At smaller redshifts, 
the quasars contribution to the UVB leads to Tjgm — 4 x lO'^K (Theuns et 
al., 2002). 

The cumulative UV background flux, J^UajZo), observed at the frequency 
and redshift Zo, (in units of 10"^^ erg cm~^ s~^ sr~^ ) due to photons emitted 
from redshifts between Zo and an effective emission screen located at Zgc > Zo, 
is the solution of the cosmological radiative transfer equation (Peebles, 1993; 
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Haiman, Rees & Loeb, 1997): 

where: jiy^i z) is the comoving emission coefficient (in units of 10~^^ erg cm~^ 
g-i computed at emission redshift z and photon frequency = Vo{^ ^ 

2;)/(l + Zo), re//(fo, -^o) -2^) is the effective optical depth at the frequency 
due to the absorption of the residual gas in the IGM between Zo and z and 
{dt/dz)~^ = -Hq{1 + z)^JVL^{l + zf + is the line element in our ACDM 
cosmological model. 

Above the hydrogen ionization threshold (i/^^ = 13.6 eV) the UV radiation 
background is processed due to the absorption of residual gas in the IGM 
dominated by neutral hydrogen and helium. At these frequencies the effective 
optical depth is given by (see e.g. Haiman, Rees & Loeb, 1997): 

[ dt 

Teffii^o, Zo, ^) = c y -^f^i^z, z)dz, k{u^, z) ^ cri{j^z)ni{z) , (2) 

where: i — {Hj , Hei) and cTj and are the cross section and the number 
density respectively, corresponding to species i. 

Below hydrogen ionization threshold, in the frequency range 11.18-13.6 cV 
(enclosing Lyman and Werner bands of H2), the UV background spectrum is 
processed due to the absorption into H2 lines (Ciardi, Ferrara & Abel, 2000; 
Haiman, Abel & Rees, 2000). The dissociation of H2 molecules removes pho- 
tons from the UVB and decreases the rate of dissociation process inside the 
collapsed halos. The magnitude of this process depends on the H2 fraction and 
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the redshift distribution of the sources. Haiman, Abel & Rees (2000) found 
that the intergalactic H2 has a neghgible effect on the UVB spectrum, both 
because its effective optical depth is small (<0.1) and because it is photodis- 
sociated at early stages. 

To compute the spectrum of the UV radiation background taking into account 
the radiative transfer effects as given by eq. (1), we assume that the comoving 
emission coefficient jli'z, has a power- law spectrum with redshift dependent 
normalization, j2i{z), and a free spectral index, a, of the form: 



The frequency range corresponds to photons collected at observed frequencies 
in 11.18 - 13.6 eV bands (LW), observed redshifts in the interval 5 - 50, arriving 
from emission screens uniformly distributed in space, located between Zo and 
a maximum redshift 2;^"^=200. 

Similar to the model in Gnedin and Hui (1998), j2i{z) was parameterized as 
a function of the reionization redshift, Zre, and two free parameters, A and B, 
in the form: 



For z Zre, the redshift dependent normalization, j2i(-2), has a value propor- 
tional to A and gradually drops for z > z^e in a redshift interval A2; ~ B~^. 
This transition period is motivated by numerical simulations (Gnedin and 




(3) 




(4) 
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Fig. 1. Panel a): Redshift evolution of the averaged UVB flux (units of 10~^^ erg 
cm~^ sr~^ ) in LW bands. Panel b): Redshift evolution of the averaged effective 
optical depth in LW bands assuming a constant H2 fraction of 2 x 10~^. The model 
parameters are indicated in the Table. 




Ostriker, 1997) showing that the hydrogen neutral fraction drops steadily 
as the formation rate of collapsed regions of high density increases in time. 
The ionizing intensity climbs gradually until the gas becomes highly ionized 
{Hj ^ 10""^), thus completing the reionization process. 

Each choice of the parameter vector p = (A, B, Zre, a) in eq. (4) defines an UV 
background flux model. 

Panel a) from Fig. 1 presents the evolution with redshift of the UV back- 
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ground flux averaged in LW bands obtained for different parameterizations of 
the comoving emission coefficient as indicated in tlie Table. Panel b) of the 
same figure presents the evolution with redshift of the effective optical depth 
averaged in LW bands computed by using the same procedure as Haiman, 
Abel & Rees (2000), for a constant H2 fraction of 2xl0~^ (r^ the H2 fraction 
at recombination). As at smaller redshifts the H2 fraction is substantially re- 
duced, the effective optical depth of the IGM due to absorption into H2 lines 
is even smaller. We found that by neglecting the absorption into H2 lines the 
UVB flux is affected by maximum 10% for redshifts higher than 30. 
Panel a) from Fig. 2 presents the evolution with redshift of the UVB flux 
at v^j}^ — 13.6 eV (?al Ryd). We also indicate (thick dotted hue) the crit- 
ical UVB flux below which (T^j^ < 10^ *^K) the star formation is prevented 
(Haiman, Abel & Rees, 2000). 

The photo-ionization rates, P-y., and the photo-heating rates, e-y., of any species 
i are related to the UVB flux, J(i/, z), through: 



^ J hu 'J hu 



where (Ti{v) are the photo-ionization (dissociation) cross-sections (Abel et al. 
1997) and vf^ are the ionizing threshold frequencies of species i (i/|^^=13.6eV, 
i/*^,^=24.6eV, i/*^,^^=54.4eV). 

We compute P^. and e-y. for those choices of parameters p = (A, B, Zre, ol) 
in eq. (4) that constraint the hydrogen photo-ionization rate at 2; ~ 6 to 
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Fig. 2. Left panel: Evolution with redshift of the UVB flux at i/^^ = 13.6 eV {^1 
Ryd). The critical UVB flux below which {Tyir < 10^-^K) the star formation is 
prevented (thick dotted line) is from Haiman, Abel Sz Rees (2000). Right panel: 
Redshift evolution of the hydrogen photo-ionization rate for the UVB models pre- 
sented in the right panel. The observational upper limits based on Lya, Ly/3 and 
Ly7 Gunn- Peterson troughs at 2; = 6.05 are from Fan et al. (2002). We indicate by 
j^thm hydrogen photo-ionization rate obtained by neglecting the optical depth 
effects for model M5. The model parameters are indicated in the Table. 




r^^^ ~ 8x10"^"^ as indicated by the observational upper limits based on 
Lya, Ly/3 and Ly7 Gunn- Peterson troughs (Fan et al. 2002). 
Panel b) from Fig. 2 shows the evolution with redshift of hydrogen photo- 
ionization rates corresponding to the UVB models presented in panel a). We 
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Fig. 3. Redshift evolution of the mean excess energy of Hj , Hei and Hen ionizing 
photons for simulation Mi (solid line) and M2 (dashed line). 
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also indicate by Mg redshift evolution of the hydrogen photo-ionization rate 
obtained by neglecting the optical depth effects for simulation M5. The com- 
parison between two cases shows that the optical depth effects results in a 
gradually increase of photo-ionization rate over a larger redshift interval. 
The differences between various UVB models can be understood in terms of 
energy input of the ionizing photons before and during the reionization. In 
Fig. 3 we show the redshift evolution of the mean excess energy of ionizing 
photons, Epf^ — e-y./r^., for models Mi and M4. Before and during reioniza- 
tion the mean excess energy is higher for model M4. For this case the redshift 
interval, Az ~ B"^, in which the reionization is completed is larger (see also 
Fig. 2). For both models the mean excess energy is dominated by the input 
energy from Hen ionizing photons. 

We find that, at 2; ~ 6 the mean excess energy from Hen ionizing photons is 
~ 0.42i/)^g for model Mi and ~ 0.5z/|^g for model M4. This corresponds to 
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a Hen photo-heating rate about 1.2 times larger for model M4 at 2; 6. The 
increased Hen photo-heating rate results in an increase of the IGM tempera- 
ture before and during the reionization, as it has been demonstrated by Abel 
and Haehnelt (2000). 

3 The ionization fraction in the presence of UVB flux 

We modify the Smooth Particle Hydrodynamics (SPH) code HYDRA(v4.2) ^ 
(Couchman, Thomas and Pearce, 1995; Theuns et. al., 1998) to study the red- 
shift evolution of the ionization fraction in the presence of UV radiation back- 
ground taking into account the radiative transfer effects, as given by eq. (1). 

Cosmological hydrodynamic simulations have proven successful in studying 
the IGM reionization at sub-galactic scales (see e.g. Hernquist et al. 1996; 
Miralda-Escude et al. 1996; Ciardi et al. 2000; Ciardi, Stoehr & White, 2002). 
It is useful to note the importance of the resolution and the volume of simu- 
lations aimed the study of reionization process. The resolution must be high 
enough to follow the formation and evolution of the objects responsabile for 
producing the bulk of the ionization radiation. At the same time, a large 
simulation volume is required in order to have a region with representative 
properties and to avoid the bias due to the variance on small scales. 
^ http: / /hydra.mcmaster.ca/hydra 
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Fig. 4. The mutual dependence on the simulation box length of the mass per gas 
particle, M^, (solid line) and of the rms linear overdensity of a sphere with the mass 
equal to that of the simulation box at present time, ctq, (dashed line). 
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Ciardi et al. (2000) derived the minimum mass of the objects that contribute 
substantially to the reionization process. They show that at 2; > 15 the main 
contribution to Hh filling factor comes from small mass objects, M ~ 10 ''M©, 
while at lower redshifts, when the formation of such objects is suppressed 
by the feedback mechanisms, the main contribution comes from objects with 
masses M > IO^Mq. 

In Fig. 4 we show the mutual dependence on the simulation box length of the 
mass per gas particle, M;,, and of the rms linear overdensity of a sphere with 
the mass equal to that of the simulation box at present time, ctq. One should 
note that once the value of ctq for a given run is selected, the redshift values of 
simulation outpus with a lower than cxo are uniquely specified [a ~ cxo/ (I + 2;)]. 
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In this paper we present cosmological hydrodynamical simulations obtained 
for 128'^ dark matter particles and the same number of gas particles in boxes 
with comoving lengths of 16h~^ Mpc and a total mass per cell of M = 
6.1 X lO^h-^MQ {Mb = 1.02 X lO^h'^MQ for gas particle and Mcdm = 5.1 x 
10^h~^MQ for cold dark matter particles). 

This choice allows us to study large enough simulated regions and to resolve 
objects with masses M ~ 5 x IO^Mq, responsable for photon production at 
reionization redshifts ~15, studied in this paper. At the same time, this choice 
ensures that at the end of the simulation the volume is large enough to provide 
a proper average (e.g. cr ~ 0.07 at 2; = 6). 

The initial particle positions and velocities were set up at Zi — 50 by using 
the transfer functions calculated with CMBFAST code (Seljak & Zaldarriaga, 
1996) for our ACDM cosmological model. 

As at early times the Compton cooling is very efficient due to the coupling of 
the free electrons (left from recombination) to the CMB photons, we assume 
an initial IGM temperature T ~ Tcmb = 2.726 x (1 + 2) ~ 139K (Theuns et 
al. 1998). 

Our cosmological simulations include the physical processes relevant for the 
primordial gas dynamics in the photo-ionized IGM (Anninos et al. 1997; Abel 
et al. 1997; Theuns et al. 1988 and references therein). The baryonic gas is 
composed by neutral (Hj , Hei ) and ionic (Hh , Hen , Hem ) species, with a pri- 
mordial composition of 76% hydrogen and 24% helium by mass. The fractional 
densities of the various gas species are denoted with the standard nomenclature 
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by normalizing their comoving number densities to the present hydrogen num- 
ber density n^, [ e.g. Hj =nHi/nH, where nH=1.88xlO"''(fib/i^/0.022) cm"^]. 
The coohng model includes coUisional ionization, recombination, coUisional 
excitation, bremsstrahlung, inverse Compton cooling on the CMB, and cos- 
mological expansion cooling. Also, the relevant heating mechanisms, photo- 
ionization (dissociation) and photo-heating of hydrogen and helium are in- 
cluded. 

HYDRA(v4.2) uses functional forms of coohng and recombination rates based 
on the atomic physics coefficients collected by Cen (1992) for use in cosmo- 
logical hydrodynamic simulations (see Appendix B from Theuns et al. 1988). 
We check out them satisfactorily against those quoted using updated atomic 
data by Abel et al. (1997). 

The simulations are integrated in time up to a final epoch, zj, where the bary- 
onic gas reach its termal equilibrium determined by the balance between the 
heating, Hp^, and cooling, C, rates defined as: 

11 

Hph = {Hie^HI + H^I^jHel + H^II^lHell)/ ' = ^ ^(7") + Cad{T) , 

i=l 

where: e^^^ , e^^ej ^iHejj photo- heating rates of Hi , Hei and Hen respectively 

and Ci are the cooling rates whose functional dependence on temperature are 
given in Appendix Bl from Theuns et al.(1998). 

In the above equation Cad — ^kBTH{z)p/iJ,mp is the adiabatic cooling rate 
due to the cosmological expansion, H{z) is the Hubble expansion rate, ks is 
the Boltzmann constant, n is the mean molecular weight and nip is the proton 
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Fig. 5. Bottom: Redshift evolution of the heating, Hph, and cooling, C, rates (units 
of erg s~^ cm~^ ) for simulation M^'^^ (thin solid and dashed lines) and M5 (thick 
solid and dashed lines). Top: Coupled evolution of the IGM temperature T and 
ionization fraction Xe obtained for simulation Mg'*'" (thin solid and dashed lines) 
and M5 (thick solid and dashed lines). 




mass. 

The code chooses the output time steps determined by the maximum instanta- 
neous values of particle velocities and accelerations. An optimal low-order inte- 
gration scheme is used for advancing particle positions and velocities (Couch- 
man, Thomas, & Pearce, 1995). 

At each simulation output time step, for each gas particle i, we read out the 
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density, Aj = Pi/pb (in units of the mean baryon density, pb) and temperature, 
Tj, (converted from the internal energy). Then, we compute the abundances 
of all neutral and ionic species and the ionization fraction defined as: 

^^^^^ e{z) + Y{i{z)+mi{z) ' 



(6) 



where e=Hii +Hen +Heni is the total free electron fraction (e^ng/n^). With 
this definition — > 1 for a fully ionized gas and x — ^ for a fully neutral gas. 
Fig. 5 presents the redshift evolution of the heating and cooling rates (bottom 
panel) and the coupled evolution of the IGM temperature and ionization frac- 
tion (top panel) for simulation M5. In order to stress the importance of the 
radiative transfer effects, we present these dependences obtained by neglecting 
the optcal depth effects (thin case) and by considering the optical depth effects 
(thick case) . The increase of photo-ionization and photo- heating rates due to 
the the optical depth effects results in an increase of the IGM temperature, 
leaving distinct features on the temporal evolution of the ionization fraction. 
In Fig. 6 we present a projection of the ionization fraction in log Aj - Tj plane 
at 2; = 6 for simulation M5. 

The ionization histories obtained for the UV radiation background models 
discussed before are presented in Fig. 7. The model parameters are chosen 
to investigate the effects of Thomson optical depth in the range Te5=0.05 - 
0.1 and to study the possibility to distinguish between models with the same 

values of t^s- Models M4 and M5 are chosen to have the same r^s and the same 
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Fig. 6. A projection of the ionization fraction in logAj - Tj plane at z = 6 from 
simulation M5. 
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value of the reionization redshift (see also the Table) . We note that while dif- 
ferent ionization histories can be distinguished by the corresponding value of 
Tes, the same value of Tgs can be produced by different parameterizations of 
UVB flux. 



4 CMB angular power spectra in presence of the UV background 

To evaluate the effects of the radiative feedback on the CMB angular power 
spectra we modified the CMBFAST code to include the redshift evolution 
of various ionization histories presented in Fig 6. The corresponding CMB 

anisotropy temperature, Ct{1), and E-modc polarization, Ce{1), power spec- 
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Fig. 7. The evolution with redshift of the ionization fraction obtained for the UVB 
models presented in Fig. 2. The model parameters and the corresponding electron 
optical depth to the Thompson scattering, Tes, are indicated in the Table. 
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tra are presented in Fig. 8. 

One should note in Fig. 8 that while Ct{1) power spectra are almost degener- 
ated, the differences in different reionization histories produce undegenerated 
signatures on Ce{1) power spectra at low multipoles {I < 50). This can be 
explained by the fact that while polarization is projecting from the epoch of 
reionization at angular frequencies / = k{riQ — rjri) ( here k is the wave number, 
770 and r]ri are the conformal times today and at the epoch of reionization) the 
temperature is projecting from the (further) last scattering surface (Zaldar- 
riaga, 1997). 

The differences between the reionization models in comparison with the ex- 
pected sensitivity of the future Planck mission can be expressed in terms 
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Fig. 8. The CMB angular power spectra, Ct{1) and Ce{1) for the reionization his- 
tories presented in Fig. 7. The model parameters and the corresponding electron 
optical depth to the Thompson scattering, Tgs, are indicated in the Table. 




of relative difference between the power spectra Ce(/) of the E polarization 
component as (Naselsky and Chiang, 2004): 



2[CeM)-Ce,M 
Ce,S) + CeAI) ' 



where the indices i and j denote different reionization models. We compare 
the amplitude of the function Dij{l) with the expected relative error of the 
Ce{1) anisotropy power spectrum for the Planck experiment. If systematic 
and foreground effects are successfully removed, the corresponding error bar 
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is given by: 



^ = -^=^—[i+n.-'c,\i)wr^)] , (8) 

where /^^j^ ~ 0.65 is the sky coverage, w = {crp9pwHM)~^, c"p is the sensitiv- 
ity per resolution element Ofwhm x Ofwhm, Wi — exp[— Z(Z + 1)/2Z^] is the 
beam window function and Ig — \/2h\20p]^jj}^. For all the frequency channels 
of the Planck experiment Ofwhm is less then ~ 30 arcminutes and at low 
multipoles {I ^ 50) the dominant term in eq. (8) is the first term, due to the 
combined cosmic and sampling variance. 

In Fig. 9 we plot the function A,j(0 obtained for the models with the same 
value for r^g compared with the Planck sensitivity at the corresponding mul- 
tipoles. One can see that in all cases the deviation function Dij{l) lies above 
the cosmic variance region, indicating that the upcoming Planck sateUite 
will allow to distinguish among E-polarization power spectra with the same 
electron optical depth and also the same reionization redshift Zre, provided 
that they have been generated by sufficiently different reionization histories. 
We observe that, in practice, the possibility to determine the CMB angular 
power spectra at Z ^ 50 with a sensitivity limited essentially only by cosmic 
and sampling variance relies on the accuracy of the removal of systematic ef- 
fects (see, e.g., the review by Mennella et al., 2004 and reference therein) and 
foreground contamination. In particular, even for an ideal experiment (i.e. 
without systematic effects) Galactic foregrounds represent the main limita- 
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Fig. 9. Relative differences between couples of E-mode polarization angular power 
spectra corresponding to models having the same value of Tes- D12 for models Mi 
and M2 (solid line) , D34 for models M3 and M4 (dashed line) and D^e for models M5 
and Mq (dash-dotted line). The model parameters and the corresponding electron 
optical depth to the Thompson scattering, Tgs, are indicated in the Table. 




tion in the recovery of the CMB angular power spectra. The accuracy of their 
modeling and subtraction from microwave maps depends on the experiment 
sensitivity and resolution. This aspect is particularly crucial for CMB polar- 
ization anisotropy analyses, because of the low polarization degree of the CMB 
anisotropy. Baccigalupi et al. (2004) show that at the Planck sensitivity level 
a blind (i.e. not critically dependent on the assumptions on the properties of 
the maps to be recovered) and accurate (i.e., in this context, precise enough 
to avoid significant errors in the recovery of the angular power spectrum of 
the CMB component) separation between the Galactic and CMB polarized 
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signals from multifrequency maps is possible at multipoles Z ^ 100 while a 
worse sensitivity has a critical impact on the accuracy of this separation. For 
this reason a high sensitivity experiment, like Planck, is crucial to remove 
all systematic effects and separate the CMB from the foregrounds with an 
accuracy good enough to measure the CMB E-mode signal at the hmit of the 
cosmic variance. 

5 Conclusions 

In this paper we study the role of radiative effects for the reionization history 
of the universe and address the detectabihty of their features by the future 
CMB polarization measurements. 

We present cosmological hydrodynamical simulations at sub-galactic scales 
that allow to resolve objects with masses of M ~ 5 x 10^ M©, responsable for 
the photon production at reionization redshifts Zre ~ 15, ensuring in the same 
time a proper average of the simulation the volume at smaller redshifts. 
Our cosmological simulations include the radiative mechanisms relevant for 
the primordial gas dynamics: photo-ionization, photo-heating and cooling of 
the hydrogen and helium in the expanding universe. 

We compute the mean specific UVB fiux as solution to the radiative transfer 
equation for those parameterizations of the comoving emission coefficient that 
constaine the hydrogen photo-ionization rate value to the upper limit indi- 
cated by the experimental measurements. 
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We find that the increase of photo-ionization and photo-heating rate due to 
optical depth effects results in a significantly contribute to the heating of the 
IGM before and during the reionization. 

We show that the main effect of the UV radiation spectrum on the temporal 
evolution of the ionization fraction is given by the value of the reionization 
redshift, Zre, and the redshift interval, A2;, in which the reionization is com- 
pleted. As reionization proceeds, the excess energy of ionizing photons build 
up a sufficient UV radiation background that rises the IGM temperature. The 
net effect is a decrease of the global ionization fraction with cosmic time over 
a limited period. 

We evaluate the effects of the UV radiation background on the CMB angular 
power spectrum taking into account different temporal evolutions of the ion- 
ization fraction. 

We find that for reionization models with degenerated CMB temperature 
anisotropy power spectra, Ct, the radiative feedback mechanisms leave dis- 
tinct signatures on the E-mode polarization power spectrum, C^, at large 
scales (/ < 50). We show that through E-mode CMB polarization power spec- 
trum measurements, the Planck experiment will have the sensitivity to dis- 
tinguish between different reionization histories even when they imply the 
same optical depth to electron scattering and degenerated Ct power spectra. 
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Table 1 

Parametrization of the comoving emission coefficient for the UVB flux models stud- 
ied in this papers. We also indicate the corresponding electron optical depths to the 
Thomson scattering,res . 



Model 


P = (A,B,Zre,a) 




Ml 


(0.065, 11.25, 12, 1.8) 


0.05 


M2 


(0.068, 12.50, 17, 2) 


0.05 


M3 


(0.075, 13.75, 15, 2) 


0.08 


M4 


(0.039, 10.00, 15, 1) 


0.08 


Ms 


(0.057, 7.70, 18, 1.8) 


0.1 


y^thin 


(0.012, 7.70, 18, 1.8) 


0.04 


Me 


(0.050, 6.90, 17, 1.8 ) 


0.1 
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